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ABSTRACT 

Exotic cosmologies predict variability of the fundamental physical constants over the 
cosmic time. Using the VLT/UVES high resolution spectra of the quasar Q0347-3819 
and unblended electronic-vibrational-rotational lines of the H2 molecule identified at 
Zahs = 3.025 we test possible changes in the proton-to-electron mass ratio /iq = mp/rric 
over the period of 11 Gyr. We obtained a new constraint on the time-averaged 



variation rate of /xg of |/i//Zo| < 5 x 10" 



'15 



{la c.l.). The estimated la uncertainty 



interval of the |A/x/yUo| ratio of about 0.004% implies that since the time when the 
H2 spectrum was formed at Zabs — 3.025, fiQ has not changed by more than a few 
thousands of a percent. 
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1 INTRODUCTION 

Kaluza-Klein (KK) type models (Super-symmetric Grand 
Unification Theory, Super string models, etc.) unify gravity 
with other fundamental forces. These models predict varia- 
tions of the fundamental physical constants over the cosmo- 
logical evolution (for a review, see, e.g., Okun' 1991). Vari- 
ations of the coupling constants of strong and electroweak 
interactions would affect the masses of the elementary par- 
ticles in a way which depends on the adopted scenario for 
the expanding universe. The model parameters of these the- 
ories can be constrained by observations. One possibility to 
test values of the physical constants at different cosmological 
epochs is to study high resolution spectra of extragalactic 
objects (Savedoff 1956). 

Such analysis has been recently carried out by a num- 
ber of authors (see, e.g., Varshalovich, Potekhin & Ivanchik 
2000, and references cited therein), and a possible variation 
of the fine-structure constant, ao — /Tic, at a level of 
Aa/ao = (-0.72 ±0.18) x 10~^ was announced by Webb et 
al. (2001) who analyzed fine-splitting lines in quasar spectra. 

The first restriction on the variability of the proton-to- 



electron mass ratio^, /xq = mp/mo, stemming from quasar 
spectra |A/i//io| < 0.13 (la c.l.) was obtained by Pagel 
(1977) who compared the observational wavelengths of Hi 
and metals as previously proposed by Thompson (1975). 
However, the derived upper limit on ^0 depends on the as- 
sumption that all elements have the same fractional ion- 
ization ratios and trace the same volume elements along 
the respective line of sight. This assumption may not be 
true in general especially for QSO absorbers where complex 
absorption- line profiles are observed at high spectral resolu- 
tion. 

The proton-to-electron mass ratio can be estimated 
more accurately from high redshift molecular hydrogen sys- 
tems. With some modifications, such measurements were 
performed for the Zabs = 2.811 H2 system from the spec- 
trum of PKS 0528-250 by Foltz, Chaffe & Black (1988), 
Varshalovich & Levshakov (1993), Cowie & Songaila (1995), 
and by Potekhin et al. (1998) who set the most stringent 
limit of |A^//io| < 1-8 x 10"" (Icr c.l.). 

In this paper, we present a new upper limit on the varia- 
tion rate of the proton-to-electron mass ratio obtained from 
the analysis of a new H2 system found at Zabs = 3.025 to- 
ward the quasar Q0347-3819. 



* Based on public data released from UVES Commissioning at 
the VLT Kueyen telescope, ESO, Paranal, Chile 
t On leave from the loffe Physico-Technical Institute, Russian 
Academy of Sciences, St. Petersburg 

I Affiliated to Observatoire de Geneve, CH— 1290 Sauverny, 
Switzerland 



The present value of the proton-to-electron mass ratio is /^o = 
1836.1526670 (Mohr k. Taylor 2000). 
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Figure 1. H2 absorption features associated with the Zabs = 3.025 damped Lyo system toward Q0347— 3819 (renormahzed intensities 
shown by dots with Icr error bars were calculated for the mean values of the local continuum deviations AC/C listed in Table 1). The 
zero radial velocity is fixed at = 3.024895. Smooth curves are the synthetic H2 profiles found by the least-squares procedure as 
described in the text. The best (x^/degree of freedom) values and the mean S/N per resolution element are shown in the corresponding 
panels. 



2 OBSERVATIONS 

High resolution spectra of the quasar Q0347-3819 were 
obtained during UVES commissioning at the VLT 8.2m 
ESO telescope and are described in detail by D'Odorico, 
Dessauges-Zavadsky & Molaro (2001) and Levshakov et 
al. (2002, LDDM hereinafter). The spectrum resolution is 
FWHM ~ 7.0 km s"^ and 5.7 km s"^ in the UV and near-IR 
ranges, respectively. The signal-to-noise ratio per resolution 
element in the UV range is S/N ~ 10 — 40. 

We identified more than 80 H2 molecular lines in a 
damped Lya system (DLA hereinafter)]^ at Zabs = 3.025 
toward Q0347-3819 (LDDM). Some of them are not suit- 
able for further analysis due to Hi Lya forest contamina- 
tion. However, we selected 15 unblended H2 lines (shown in 

^ The damped Ly« systems are believed to originate in the inter- 
vening galaxies or proto-galaxies located at cosmological distances 
(e.g. Wolfe et al. 1995). 



Fig. 1) which provide the most accurate line center measure- 
ments to set an upper limit on possible changes of /io- 

We would like to emphasize that our analysis of the 
Zabs = 3.025 DLA has shown that all the H2 line profiles 
can be adequately described with a unique value of — 
3.024855 ± 0.0000050 which implies that no assumption on 
the variability of /iq is needed or statistically justified. Ob- 
servations show that the ratio A/j./ fio = {yi-z ~ ^J^o) /mo is zero. 
The analysis presented below has the objective to estimate 
at what accuracy Ap/yUo = 0. 



3 DATA ANALYSIS AND RESULTS 

For measurements of the absorption line centers in QSO 
spectra, there are three principal sources of statistical errors 

^ Based on the H2 laboratory wavelengths given in Abgrall & 
Roueff (1989). 



m-p/me non-variability 3 



Table 1. H2 lines at z^bs = 3.025 toward Q0347-3819 and sensitivity coefficients K. 



J 


Line 


Ag, A 


AC/C^ % 


-^obsi A 


z 


K 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


W3-0Q 


947.4218±0.0005 


0.9 ±0.2 


3813.2653 ±0.0041 


3.024887(5)'' 


0.0217427(8) 




Ll-OP 


1094.0522±0.0051'= 


0.8 ±0.1 


4403.4575 ± 0.0077 


3.02491(2) 


-0.0023282(1) 




L3-0R 


1063.4603±0.0001= 


0.8 ±0.1 


4280.3051 ± 0.0030 


3.024885(3) 


0.0112526(4) 




L3-0P 




0.9 ± 0.2 




3 n248Q4C4") 


m 02682(^4") 




L7-0R 


1013.4412±0.0020 


0.8 ±0.1 


4078.9785 ± 0.0025 


3.024898(9) 


0.03050(1) 




LIO-OP 


982.8340±0.0006 


0.8 ±0.1 


3955.8049 ± 0.0051 


3.024896(6) 


0.04054(6) 


2 


WO-OR 


1009. 0233±0. 0007 


0.8 ±0.1 


4061.2194 ±0.0061 


3.024901(7) 


-0.0050567(7) 




WO-OQ 


1010.9380±0.0001 


0.8 ±0.1 


4068.9088 ± 0.0053 


3.024885(6) 


-0.0068461(2) 




Wl-OQ 


987.9744±0.0020 


0.9 ±0.1 


3976.4943 ± 0.0049 


3.02490(1) 


0.0039207(2) 




L3-0R 


1064.9935±0.0009 


1.0 ±0.2 


4286.4755 ± 0.0098 


3.02488(1) 


0.0097740(5) 




L4-0R 


1051.4981±0.0004= 


0.8 ±0.1 


4232.1793 ± 0.0085 


3.024904(8) 


0.015220(1) 




L5-0R 


1038.6855±0.0032 


0.8 ±0.1 


4180.6048 ± 0.0070 


3.02490(1) 


0.020209(3) 


3 


L4-0R 


1053. 9770±0. 0011 


0.8 ±0.1 


4242.1419 ±0.0028 


3.024890(5) 


0.012837(2) 




L6-0R 


1028. 9832±0. 0016 


0.8 ±0.1 


4141.5758 ±0.0031 


3.024921(7) 


0.022332(7) 




L12-0R 


967.6752±0.0021 


0.8 ±0.1 


3894.7974 ±0.0039 


3.02490(1) 


0.0440(2) 



Notes: "listed values are from Abgrall et al. (1993a, 1993b); ''the local continuum deviation with la 
error; the error is estimated from the comparison between the data in Abgrall & Roueff (1989) and 
Abgrall ct al. (1993a, 1993b): la standard deviations are shown in parenthesis (last digits after a 
decimal point), e.g. 3.024887(5) means 3.024887 ± 0.000005. 



caused both by the data reduction procedure and by sta- 
tistical fluctuations in the recorded counts. First, all echelle 
data must be rebinned to constant wavelength bins in order 
to combine different spectra and hence to increase the signal- 
to-noise ratio. Such resampling changes statistical properties 
of the noise and introduces correlations between the data 
point values. The correlation coefBcient of about ±0.8 was 
estimated by LDDM for the UVES data in question. This 
means that the classical test is no longer applicable to 
the rebinned spectra, because the measurements are not in- 
dependent. The second error, usually important for lines ob- 
served in the Lya forest, arises from the uncertainty of the 
overall continuum level and its shape in the vicinities of in- 
dividual spectral features. The third error is connected with 
the finite number of photons counted. The individual inten- 
sity points within the interval of integration can fluctuate 
producing distortions in the line shape. 

There are, of course, other sources of errors in the mear 
surements of the precise wavelengths Aobs of the observed 
H2 lines. The values of Aobs may be affected by errors in 
the wavelength calibration and by the thermal shifts. The 
observed wavelengths are not significantly affected by errors 
in the correction of the wavelength scale to the heliocentric 
and vacuum values. Typically, these errors do not exceed a 
few m s~^ (e.g. Edlen 1953). The quasar spectrum is, how- 
ever, constantly shifted with temperature by 0.37 pixel/K. 
Given the approximate pixel size (0.02 A) at 4000 A, we 
can expect the shift of at most 0.3 km s~^ in case the Th 
calibration spectrum was taken a few hours after the science 
exposures and the difference in temperature was 0.5 K. The 
thermal shift should be more or less constant for all lines, so 
it affects only the assumed redshift. The values of Aobs are 
mainly affected by errors in the wavelength calibration. We 
found that the rnis of this systematic error over the total 
range covered by the blue spectrum is Usyst = ±0.0016 A. 
This corresponds at 4000 A to 0.12 km s~^ 

The continuum placement in the UV portion of the 



Q0347-3819 spectrum was determined by using 'continuum 
windows' in the Lya forest which were fitted by a low order 
polynomial. The accuracy 5c of the local continuum (which 
may deviate from the common continuum in the Lya forest 
region) was estimated during the fitting procedure described 
below. 

For a given H2 line, the number of pixels involved in 
the analysis corresponds to the number of points in the line 
profile shown in Fig. 1. To measure the line centers we used 
a method that matches the observed profiles with the syn- 
thesized ones to estimate the set of model parameters. Our 
previous analysis of the line profiles of different elements 
identified in the 2abs = 3.025 H2-bearing cloud has shown 
that both low ions and H2 have the same simple velocity 
structure of the main component - a narrow symmetrical 
core with a common broadening parameter = 2.80±0.45 
km s~^ (LDDM). Thus in our present study we also used a 
simple one-component model with four free parameters: the 
center, the width, the intensity of the absorption line and 
the local contirmum displacement Sc — AC/C (this tech- 
nique was successfully used in the study of metal lines in 
the Lya forest by Molaxo et al. 2001). The set of initial 
parameters then was adjusted until a satisfactory fit could 
be achieved. The objective function was augmented with 
the penalty function in the form [cf. Eqs. (13) and (14) in 
LDDM]: 



b - feHa 



(1) 



where fena = 2.80 km and o-t^^ = 0.45 km s~^. 

To evaluate statistical errors, Monte Carlo analysis was 
performed in the same way as described in Appendix A in 
LDDM. The results of these calculations are presented in 
columns (4) and (5) of Table 1. It is seen that for the selected 
H2 lines 6c ^ 1%. This means that the H2 lines involved in 
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• J = 1 

■ J = 2 

i J = 3 



-0,01 0.01 0.02 0,03 0,04 0,05 

Sensitivity Coefficient, K 

Figure 2. Relation between tiie rcdshift values Zi calculated for 
individual H2 features shown in Fig. 1 and their sensitivity coeffi- 
cients Ki. The solid line is the linear regression: Zi = z + niJCi — K.) 
with K = (1 + z)Afi/fio. The dotted lines representing the la de- 
viations from the slope k of the best linear regression demonstrate 
that A/i/iMo = at the level ~ 3 X 10~^. 



■ (J=1)+(J = 2)+CJ = 3) 
(J = 1)-(J = 2) 




(Am/Mo)x10= 

Figure 3. Monte Carlo simulations of the probability density 
functions (p.d.f.) of A/i/fio for two samples: (i) the H2 lines 
arising from the low rotational levels with J = 1,2,3 (the solid 
curve), and (ii) the H2 lines arising from the low rotational levels 
with J = 1, 2 (the dashed curve). The difference between p.d.f. is 
caused by gradual shift in the radial velocity as J increases. 



the analysis are not significantly disturbed by the Lya forest 
absorption. 

In this approach the central wavelengths Aobs can be 
measured with an accuracy exceeded the average pixel size. 
As expected, the weaker absorption features have lower ac- 
curacy compared with the stronger lines. The correspond- 
ing redshifts for individual H2 lines are shown in column 
(6) with the la uncertainties including osyst- In these cal- 
culations we used the II2 laboratory wavelengths Ao from 
Abgrall et al. (1993a, 1993b) which show slightly diflerent 
Ao values as compared with their previous publication (Ab- 
grall & Roueff 1989). Therefore, as a consequence the new 
values of discussed below are shifted with respect to 
those presented in LDDM. 

With the obtained z values and their standard devia- 
tions we are able to constrain possible changes of fio- Indeed, 
if the ratio Afj./fj,o = (/iz — fxo)/ jjbo is not zero, then for any 
two H2 lines with rest frame wavelengths Ai,o and Aj,o the 
ratio 

^i:i/^^l-f(/C,-/C, ) Am/mo (2) 

would deviate from unity, where A^.^ and Aj,z are the corre- 
sponding line centers measured in a quasar spectrum. 

In equation (^), the so-called sensitivity coefficients K. 
determine the sensitivity of H2 wavelengths to the variation 
of the proton-to-electron mass ratio. These coefficients have 
been calculated by Varshalovich & Levshakov (1993) and 
in a different manner by Varshalovich & Potekhin (1995). 
Both methods give results in good agreement. The most 
recent version of the procedure to calculate /C and Afi/fio 
from QSO spectra is described in detail by Potekhin et al. 
(1998). Following this procedure we calculated the sensitiv- 
ity coefficients which are listed in column (7) of Table 1. The 
accuracy of the IC values was also estimated from coefficients 



Ymn (listed in Table 1 in Potekhin et al.) using the method 
of error propagation {Ymn values were considered to be ac- 
curate to k decimal places and their rounding errors were set 
to 0.5 X 10"*°). It should be noted that although transitions 
between the excited states of the H2 molecule have higher 
wavelength-to-mass sensitivity coefficients, their accuracy is 
much lower as compared with the low- J transitions. In linear 
approximation 

Zi = z + K^fCi — K.) , (3) 

where n — (1 + z)Afi/fio with z and K, being the mean 
redshift and the mean sensitivity coefficient, respectively. 

The linear regression in the form ^ was firstly calcu- 
lated for the complete sample of the H2 lines (Table 1) where 
transitions from the J = 1,2 and 3 levels are combined. 
The obtained result (A^/^o).j=i+2-t-3 = (5.0 ± 3.2) x 10"^ 
(which is consistent with the value independently found by 
Ivanchik et al. 2002) shows a "possible variation" of fio at 
the 1.6(7 level. However, if we consider II2 transitions from 
individual J levels, then the weighted mean redshifts reveal 
a gradual shift in the radial velocity for features arising from 
progressively higher rotational levels : z^"^^ = 3.024890(2), 
zi^=2 = 3.024895(3), and z^=^ = 3.024904(4)[|. The H2 lines 
with changing profiles and small velocity shifts as J increases 
were also observed in our Galaxy in the direction of ^ Ori A 
by Jenkins & Peimbert (1997) who suggested that these H2 
lines are formed in different zones of a postshock gas. 

Albeit the Zh^^ and Zh^^^ values are consistent within 
1(7 intervals, the difference between z^"^^ and Zh2^^ is essen- 
tial and equals 1.0 ±0.3 km s"'^. If we now exclude from the 
regression analysis levels with J — 3, then (A////io) j=i-i-2 = 

^ The numbers shown in parenthesis are the la standard devia- 
tions in last digits after a decimal point. 
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(2.1 ±3.6) X 10^''. (The use of samples with H2 lines arising 
from the same rotational levels would be more reasonable 
to estimate A/i//io, but in our case the sample size is rather 
small and we have to combine the J = 1 and J = 2 levels 
to increase accuracy). This hncar regression is shown by the 
solid line in Fig. 2 while two dashed lines correspond to the 
la deviations of the slope parameter k. 

We also calculated the probability density functions of 
A/i/^o for two samples of the H2 lines (with J = 1,2,3, and 
J — 1,2) using statistical Monte Carlo simulations which 
suggest that the errors az and aic are normally distributed 
around the mean values of z and K. with the dispersions 
equal to their probable errors listed in Tabic 1. The result 
is presented in Fig. 3. It is clearly seen that small changes 
in the radial velocity with increasing J cause the difference 
between these two probability density functions and thus 
may mimic a shift in A/x//xo- 

From these calculations we find |A/i//i()l < 5.7 x 10~^ 
(la) which is about throe times stronger as compared with 
the value estimated by Potckhin et al. (1998). 

Both constraints on the variation of the m-p/rtie ra- 
tio are in good agreement with the limit on the variabil- 
ity of the product of the fine-structure constant, nuclear 
g factor of the proton and the masses of the electron and 
proton Aln (a^ (jp mc/rrip) = (1.2 ± 1.8) x 10~* which was 
set by de Bruyn, O'Dea k. Baum (1996) from the mea^ 
surements of the redshifts of the Hi 21 cm line and the 
optical resonance absorption lines observed at Zabs = 3.38 
in the DLA toward Q0201-I-113. The same order of mag- 
nitude restrictions to Aln (a^ me/mp) were found from 
other DLAs detected at 2abs = 0.524 (Q0235-I-164), 0.692 
(3C286), 1.944 (Q1157-h014), and 2.038 (Q0458-020) (see 
Table 4 in de Bruyn et al. 1996), and even stronger limit 
of (0.7 ± 1.1) X lO"'^ at 2abs= 1.776 (Q1331-I-170) was 
set by Cowie & Songaila (1995). Comparison of Hi 21 cm 
and molecular absorption (CO, ^^CO, C^*0, CS, HCO+, 
and HCN) also yields very tight constrains on the ratio 
|A(a2 Pp)/(q2 gp)o| < 0.7 X lO"-'"' at z^bs = 0.2467 and 0.6847 
toward Q1413-I-135 and Q0218-I-357, respectively (Murphy 
et al. 2001). 

In our calculations the la confidence interval to Afi/fio 

was set as 

-1.5 X 10"^ < Afj,/fj,o < 5.7 X 10"^ . 

For a cosmological model with Om = 0.3, Ha = 0.7, and 
Ho = 72 km Mpc~^, the look-back time for z^ha = 3.025 
is 11.2 Gyr [see, e.g., equation (16) in Carroll, Press, & 
Turner 1992]. This leads to the restriction 

\fl/^io\ < 5 X lO"'"^ yr"^ 

on the variation rate of po ■ 



Since the functional dependence of the masses of proton 
and electron on the fine-structure constant is unknown, we 
are not able to compare directly our result with changes in 
a at the level of ~ 0.7 x 10"^ found by Webb et al. (2001). 
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4 CONCLUSIONS 

We have obtained a new constraint on the variation rate of 
the proton-to-electron mass ratio of A/i//io = (2.1 ± 3.6) x 
10"^ at «abs= 3.025 toward Q0347-3819 {At ~ 11 Gyr). 
The accuracy is a factor of 3 higher as compared with the 
measurements in another H2-bearing cloud at Zabs = 2.811 
toward Q0528-250 (Potekhin et al. 1998). Both measure- 
ments show no statistically significant changes of fio on space 
and time coordinates. 



